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Abstract 

We introduce a model of a lossy second-harmonic-generating (x^^-*) cavity 
externally pumped at the third harmonic, which gives rise to driving terms 
of a new type, corresponding to a cross-parametric gain. The equation for 
the fundamental-frequency (FF) wave may also contain a quadratic self-driving 
term, which is generated by the cubic nonlinearity of the medium. Unlike pre- 
viously studied phase-matched models of x^^^^ cavities driven at the second har- 
monic (SH) or at FF, the present one admits an exact analytical solution for the 
soliton, at a special value of the gain parameter. Two families of solitons are 
found in a numerical form, and their stability area is identified through numer- 
ical computation of the perturbation eigenvalues (stability of the zero solution, 
which is a necessary condition for the soliton's stability, is investigated in an 
analytical form). One family is a continuation of the special analytical solution. 
At given values of parameters, one soliton is stable and the other one is not; 
they swap their stability at a critical value of the mismatch parameter. The 
stability of the solitons is also verified in direct simulations, which demonstrate 
that the unstable pulse rearranges itself into the stable one, or into a delocalized 
state, or decays to zero. A soliton which was given an initial boost C starts to 
move but quickly comes to a halt, if the boost is smaller than a critical value 
Ccr- If C > Ccr, the boost destroys the soliton (sometimes, through splitting 
into two secondary pulses). Interactions between initially separated solitons are 
investigated too. It is concluded that stable solitons always merge into a single 
one. In the system with weak loss, it appears in a vibrating form, slowly re- 
laxing to the static shape. With stronger loss, the final soliton emerges in the 
stationary form. 

PACS numbers: 42.65.Yj; 42.65. Tg; 05.45.Yv 
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1 Introduction 



In the vast family of optical solitons, an important niche is occupied by solitary 
waves in cavities, supported by the quadratic (x^^-*) nonlinearity of the degen- 
erate optical-parametric-oscillator (OPO) type 0111 ■ The intrinsic loss in the 
cavity should be compensated by an external pump field E, which in most cases 
is supplied at the second harmonic (SH) T. Such an arrangement is frequently 
referred to as downconversion, and is described by the model in which either 
the evolution equation for the SH field v explicitly contains a constant driving 
term ^ E, or, equivalently, the equation for the fundamental- frequency 

(FF) wave u includes a parametric-gain term ^ Eu* , where the asterisk stands 
for the complex conjugation |S1[71|H]- Alternatively, the cavity may be exter- 
nally pumped at the FF, which is referred to as upconversion, the respective 
model including a constant driving term in the FF equation For both cases, 
families of one-dimensional solitons and their stability have been investigated in 
detail, see Refs. 01-201 and references therein. The stability of cavity solitons 
was tested in a direct experiment using the photorefractive nonlinearity 

It should be mentioned that, as the cavity models are dissipative ones, the 
solitary pulses found in these models are not solitons in the rigorous sense. 
Nevertheless, this term is broadly applied to them, therefore we use it in this 
paper too. 

The analysis of solitons in OPO models was extended in various directions, 
including the study of moving solitons interactions between them jl3) . 

nondegenerate second-harmonic-generation (SHG) which involves two FF waves 
with orthogonal polarizations jJH, general three- wave interactions |15j, the use 
of the quasi-phase- matching technique two-dimensional solitons (see, e.g., 
Ref. 23), etc. Besides their significance as the subject of fundamental research, 
OPO cavity solitons also have a potential for the design of rewritable multi-pixel 
optical- memory patterns 

In this work, we aim to propose and analyze another possibility to drive dis- 
sipative cavities with the SHG nonlinearity, namely, through a phase-matched 
third-harmonic (TH) pump wave, w. Obviously, the corresponding driving 
terms in the equations for the FF and SH fields u and v, which are generated 
by the parametric interaction of these fields with the TH pump, are propor- 
tional, respectively, to wv* and wu* . Thus, this model includes a new feature, 
the cross-parametric gain, in the system of coupled FF and SH waves, and a 
problem of straightforward interest is to study solitons that can be supported 
by this type of the gain in the lossy SHG setting, especially as concerns stability 
of the solitons. Besides that, we will also take into regard a possibility of an 
additional quadratically nonlinear parametric self-driving term in the equation 
for the FF field u, in the form of w (u*)'^ , which may be induced by the same 
TH pump through the x^^^ (cubic) nonlinearity. In this work, we consider only 
bright solitons; dark solitons, and patterns in the form of domain walls (see, 
e.g., Ref. 19J may also be possible in the model including the TH drive. 

As concerns the mutual phase matching between the FF, SH, and TH fields, 
which is implied in the model, it was demonstrated, in another context, that 
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matching of this type may take place in the so-called multi-step x^^^ systems f2U\ . 
However, the subject of Ref. [50] was the corresponding three-wave solitons, 
rather than pumping of the FF and SH fields through the TH wave. 

The paper is organized as follows. In Section II, we give the formulation 
of the model, and find particular exact solutions for the solitons in an ana- 
lytical form. In the same section, we also investigate stability conditions for 
the zero solution, which is a necessary prerequisite for the stability of solitons. 
In Section III, we find two families of general soliton solutions in a numerical 
form. One family is a direct continuation of the exact solution, while the other 
one is different. As well as the analytically found soliton, they always have 
a single-humped shape, but, unlike the constant-phase exact solutions, they 
feature intrinsic chirp. Stability regions for the solitons are identified, in the 
system's parameter space, through computation of the corresponding eigenval- 
ues for small perturbations. One of the solitons is stable, and the other one is 
not; they swap the stability via a bifurcation at a critical value of the mismatch 
parameter. The shape of the stability regions is quite nontrivial; in the model 
with weak losses, the stable solitons may have complex eigenvalues, correspond- 
ing to weakly damped intrinsic oscillatory modes. The stability is also verified 
in direct simulations, with a conclusion that the unstable soliton rearranges into 
the stable one, or into a delocalized spatially periodic state, or, sometimes, it 
decay to zero. Also in Section III, we show that attempts to produce moving 
solitons fail: if the soliton is initially boosted by lending it a speed C, it quickly 
stops, provided that C is smaller than a critical value Ca', the boost with C > 
Ccr destroys the soliton. Section IV deals with interactions between two stable 
solitons, initially placed at some distance. They merge into a single soliton, 
which emerges with slowly fading intrinsic vibrations in the weakly dissipative 
model, or immediately in the stationary form if the loss parameter is larger. 
Section V discusses possible extensions of the work and concludes the paper. 

2 The model and exact results 

Equations which describe the degenerate x^^^ interaction between the FF and 
SH fields u{x, t) and v{x, t) in a one-dimensional lossy cavity in the presence of an 
additional pump TH wave wq, whose depletion is negligible, are straightforward 
to derive: 

{qi — iai) u + iwQV* + ix^^^'wo (u*)^ , (1) 

2 {q2 ~ ia2) V + iwQU* . (2) 

Here, loq is the fundamental carrier frequency, the coefficient amenable for the 
three- wave coupling between the FF, SH and TH fields is absorbed into wq, qi 
and q2 are real detuning coefficients at the FF and SH, and ai and a2 are the 
loss coefficients for the same fields. In this paper, we focus on the most natural 
case, ai = a2, but the analysis has demonstrated that the results do nor differ 
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in any noticeable aspect for oi\ ^012- The last term in Eq. ^ takes into regard 
the above-mentioned possibility of the nonlinear parametric self-driving of the 
FF field, which can directly couple to the TH pump wave in the presence of 
the cubic nonlinearity. The factors of i in front of the cross-parametric-driving 
terms in Eqs. ^ and Q, while is assumed real, can always be fixed, defining 
phase shifts between the corresponding waves. 

After obvious normalizations (in particular, the SH field is rescaled with 
u — *■ \/2u, so as to keep the coefficients in front of the cross-driving terms equal 
in the two equations), Eqs. (Q) and l(2Jl can be cast in the following form: 

(1 — iai) u + iaQV* + i(3 (u*)^ , (3) 

{2q — iai) V + iaou* , (4) 

where ao and /? are the effective pumping coefficients, both proportional to wq, 
the mismatch coefficient in Eq. ^ is normalized to be 1 (a different variant of 
the model is obtained by fixing the latter coefhcicnt to be —1, but the existence 
of bright solitons is not expected in that case). The notation for the fields u 
and V and the coefficient ai was not altered, although they were rescaled (recall 
we assume a2 = 0:1). By means of a phase shift of u and v, ag may always be 
made real, which we assume below, while /3 is, generally speaking, a complex 
coefficient. 

Obviously, Eqs. Q and Q cannot give rise to stable localized solutions 
unless the zero solution, w = u = 0, is stable. Linearizing the equations, an 
elementary calculation yields the following stability condition for the zero solu- 
tion: 

al < al + {l + 2qf /4, if 1 + 2? > 0; 

al < aj, if 1 + 2(7 < 0. (5) 

The consideration of the linearized version of Eqs. Q and makes it also 
possible to predict the asymptotic form of the exponentially decaying tails of 
the soliton solution at |a;| ^00, 

{u,v*)-^exp{-X\x\), (6) 

where a (generally) complex constant A (it is defined so that its real part is 
positive) is to be found from the equation 

QA^-l + zaij QA2-2g-zai^ =ag . (7) 

The complex structure of A implies that the soliton must have a nontrivial 
intrinsic phase structure, which will be studied in detail below. 

Particular exact soliton solutions to Eqs. ||2J) and Q can be sought for as 

uoix) = Ae'-^sech^ (kx) , ^0(2:) = B e^^'^sech^ (kx) , (8) 



iut + -Uxx + u V = 



IVt 
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where k, (f> and A, B are real constants. This ansatz may produce solutions in 
the case of /3 = [no self-driving quadratic term in Eq. Q]. Then, the soliton's 
parameters are obtained by direct substitution of the expressions © into Eqs. 
© and gl): 

K^ = l- 2q, (9) 

B ^ 3K^ (10) 

cos {3(f>) — ai / , sin (30) = (1 — 4^) /ao . (11) 

Due to the identity sin^ (30) -|-cos^ (30) = 1, Eqs. give rise to an additional 
constraint on the parameters of the model which is necessary for the existence 
of the exact solution in the above form, 

«o = (ao)exact = \/ + (l ^ Aqf . (12) 

Equation (|12(l determines the value of the cross-parametric-gain which is 
necessary to support the exact soliton solution. Thus, the two conditions, /3 ~ 
and H12|l , must be imposed on the parameters of Eqs. © and (0J) to provide for 
the existence of the analytical solution for the soliton in the simple form ^ (in 
fact, the condition ai =0.2-, which was adopted above, is also necessary for the 
existence of the exact soliton in the present form). 

It is relevant to notice that, comparing Eq. with the asymptotic waveform 
(O, one can identify, for the present solution, A = 2k. Then, taking into regard 
Eq. ifT^ . it is easy to verify that this A indeed satisfies Eq. {T)). 

A necessary stability criterion for the exact soliton solution can be obtained, 
inserting the relation H12|l into the stability condition for the zero solution. 
After a simple algebra, it takes the form 

1/10 <q< 1/2. (13) 

A typical example of the exact soliton given by the above expressions is 
shown in Fig. ^ Below, numerical results will be produced for ao ^ (Q^o)cxact 
and /3 ^ 0, and they will be compared to the shape shown in Fig. ^ 

It is relevant to mention that the previously considered phase-matched OPO 
models with the parametric gain of the downconversion type (provided by the 
SH pump) did not produced exact solutions similar to the present one. Exact 
solutions were only obtained in the case of large detuning In that limit, 

the SH field can be eliminated, and the remaining FF equation amounts to 
a parametrically driven damped cubic nonlinear Schrodinger equation, which 
has a pair of well-known exact solitary-pulse solutions (see, e.g., Ref. |22|'). 
On the other hand, nongeneric exact analytical solutions for solitons (which 
explicitly contain the chirp) were found in a model of a two-core system 
of the Ginzburg-Landau type, in which intrinsic gain was set in the nonlinear 
core, and a linearly coupled additional lossy core played a stabilizing role |23| . 
However, the gain in that system was not of the parametric type. 
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(a) (b) 

Figure 1: A typical shape of the soliton's FF (a) and SH components (b), as 
given by the exact analytical expressions Hl()|l . and 1)11(1 for /3 = 0, 
ai = 0.1, and q = 0.2; accordingly, Eq. ifT^ yields, in this case, (Q;o)oxact ~ 
O.lVSw 0.224. 



3 The soliton family: numerical results 

The analytical solutions were found above in the special case only, and even in 
that case, full stability analysis requires the use of numerical methods. In this 
section, we aim to construct a general family of soliton solutions in a numerical 
form, and then study their stability. A possibility of the existence of moving 
solitons will be considered too. 

3.1 Stationary solitons 

A family of soliton solutions to the stationary version of Eqs. © and with 
Ut = Vt = 0, was constructed by means of a continuation procedure (based on 
the Newton's numerical method), starting with the exact solutions in the form 
given by Eqs. lO - l(TT|l . which are valid in the case of /3 = and ao = (o^o) exact' 
and then gradually varying both ag and (3. As a result, it was concluded that the 
shape of the soliton in both the FF and SH components, \u{x)\ and |?^(a;)|, does 
not not vary much in comparison with that of the exact solution, while a new 
feature is an intrinsic phase structure of the soliton's wave field, characterized by 
nonzero chirps, 4>" and ■0", where (f'i.^) V'l^^) are phases of the fields u{x) and 
v{x). These features are illustrated by Figs. ISlandEl which display two generic 
subfamilies of the numerically found stationary soliton solutions, obtained by 
varying the self-driving coefficient (3 at fixed values of the cross-parametric-drive 
coefficient ao (in these figures, only real values of /? are included; complex (3, 
which do not produce effects drastically different from those found for real /?, will 
be briefly considered below). Figure |2] corresponds to ag = (o^o) exact' which is 
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the value that gives rise to the exact solution (for (3 = 0), see Eq. ifT^ . and Fig. 
121 presents the situation at a different value of ao, namely, ag = 0.9 (ao)exacf 
In these figures, the amplitude profiles are shown through their differences from 
those corresponding to the exact soliton solution (jSJ - taken at /? = and 
ao = (Q^o)exact' bccausc full profiles are too close to each other. 

At other values of ao, the soliton solutions are quite similar to those displayed 
in Figs. [5]and|21 We stress that, while the amplitude-profile differences displayed 
in these figures may have a relatively complex form, the full profiles of 
and \v{x)\ always keep the simple single-humped shape. 

3.2 Stability of the stationary solitons 

The most fundamental approach to the investigation of the stability of station- 
ary solitons is based on computation of stability eigenvalues within the frame- 
work of the system of linearized equations for infinitesimal perturbations. We 
have performed this computation through the corresponding Jacobian matrix of 
the linearized system. The results can be conveniently summarized in the form 
of maps showing stable and unstable regions in the system's parameter planes. 
First, in Fig. 21 we display the stability map for the exact soliton solutions © - 
((TT)l in the plane (ai, g), with (3 = and ao = (ao)exact' ^'i- 03, s.nd also 
for the solitons found numerically at the same values of the parameters, except 
for the self-driving FF coefficient, (3 = 0.1. Note that the stability region for the 
exact solitons is located inside the stripe corresponding to the necessary sta- 
bility condition (|13|l . being actually much narrower than it (which means that 
there are strong nontrivial stability conditions for the soliton proper, which do 
not amount to the simple criterion that guarantees the stability of the zero 
background). 

In Fig. 0](and similarly in Fig. El see below) we distinguish between stability 
sub-regions in which all the eigenvalues are real, and those where complex ones 
are found. Accordingly, a perturbation applied to the stable soliton in the latter 
sub-region excites a damped intrinsic oscillatory mode (this will be observed 
below, in the case of interactions between the solitons). Quite naturally, the 
complex eigenvalues are found in the case when the loss parameter ai is small 
enough. 

A similar stability map is shown in Fig. [SI for the same case which was se- 
lected above for Fig. 13 i.e., ao = 0.9 (ao)(,xacf Additionally, for both cases ao = 
("o)oxact and ao = 0.9 (ao)^^^^' as well as for another one, ao = 1.1 (ao)oxact> 
Fig. shows the area of the stability region as a function of (3, between the 
two values chosen for the display in Figs. 01 and O (3 = and (3 = 0.1. Note 
that two plots in Fig. pertain to exactly the same soliton subfamilies which 
were included in Figs. (21 and |31 A noticeable observation is that, depending 
on the value of ao, the stability area may both decrease and increase with 
As concerns the additional case of ao = 1-1 (o^o) exact' included in Fig. El the 
stability map for it is not very different from that for ao = (o^o) exact' ^^S- 
^ therefore it is not displayed here separately. 

We have also investigated the case when the self-driving coefficient (3 in the 
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Figure 2: The differences, \u{x)\ — \uo{x)\ and \v{x)\ — |'yo(a;)|, between the 
ampHtude distributions in the FF (a) and SH (b) components of the numerically 
found generic solitons and the exact analytical solution. The panels (b) and (d) 
display the phase chirp, 0" and "0", in the FF and SH components of the generic 
soliton solutions (the chirp distributions are shown on the half-axis, in view of 
the soliton's symmetry, and only in the region of x where the soliton is actually 
located). The fixed parameters are q — 0.2, ai = 0.1, and, for this figure, the 
corresponding value (ao)oxact — 0.224, as per Eq. l(T^ . is chosen. Note that the 
full shapes, \u{x)\ and |w(a;)|, of all the solitons in both components are strictly 
single-humped ones. 
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Figure 4: The stability maps for the family of the exact soliton solutions lO 
- (fTT|l . corresponding to /3 = and ao — («o)cxact' ^^1- C21), and for 
the solitons found numerically for the same ao but [3 = 0.1. The stability is 
identified from the computation of the stability eigenvalues (and verified by 
direct simulations). The shading, bordered by the dashed lines, marks stability 
sub-regions in which complex eigenvalues were found, while, in the unshaded 
areas, all the eigenvalues are real. 
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Figure 5: The same as in Fig. 0] but for ao = 0.9 (ao)^ 
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Figure 6: The area of the stabihty region vs. /3, for three different values of ag- 



FF equation © is complex. In this case, the stability maps are not drastically 
different from the ones displayed above for real /? (therefore we do not show them 
here), although the area of the stability region gets somewhat smaller. The area 
is shown, as a function of |/3|, for two different cases, /3 = (1 + «) \(}\/^/2 and 
/3 = z|/?|, in Fig.EI 
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Figure 7: The same as in Fig. Elfor ao = (Q!o)oxact ^^"^ complex (3. 



The conclusions concerning the stability of the solitons, which were drawn 
above on the basis of the stability eigenvalues, where also verified against direct 
simulations of the full system of Eqs. l|2Jl and performed at a sufficiently 
dense grid of points covering the predicted stability and instability regions in the 
system's phase space. As a result, it has been concluded that all the solitons. 
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which were predicted to be stable, are stable indeed in the direct evolution. 
The solitons which are expected to be unstable decay to zero under the action 
of small perturbations or, if the perturbation is stronger, they may rearrange 
themselves into solitons of a new kind, as described below. 

3.3 A bifurcation to the second type of stable solitons 

One may observe that the stability areas displayed in Figs. 01and[Slare located 
at q < 0.25. In particular, although we have no analytical proof of the instability 
of the exact soliton (jSJ at 5 > 0.25, we notice that the value of the mismatch 
parameter q = = 0.25 is a special one for the exact solution, as Eqs. Hll|l 
show that the phase (j) vanishes precisely at q — Qc- 

A stable soliton of a different type can be found at g > 0.25. This solution, 
which we will call a type-II soliton, to distinguish it from the one considered 
above, that we will refer to as type-I, cannot be obtained by the continuation of 
the exact solution (jSJ. This makes finding the type-II soliton directly from the 
numerical solution of the stationary version of Eqs. ^ and (|2Jl problematic, as 
a good initial guess is not available. Nevertheless, it was found, in the region 
q > 0.25, in the following way: one can take the numerically exact unstable 
type-I soliton, and add an arbitrary perturbation to it. A small perturbation 
initiates decay of the unstable soliton to zero. However, if the perturbation is 
sufficiently large, the outcome may be spontaneous rearrangement of the pulse 
into a new stable one, which we identify as the type-II soliton. An example of 
the unstable type-I soliton (alias separatrix soliton, see below), and its stable 
type-II counterpart, which is generated from it by a perturbation, is displayed 
in Fig. 

The stability map for the type-II solitons was identified, as it was done 
above for their type-I counterparts in the region q < 0.25, both by means of 
the computation of the eigenvalues from the linearized equations, and by means 
of direct simulations. The resulting map is displayed in Fig. |51 It shows not 
only the region where the type-II soliton is stable, but also indicates if the 
perturbation-induced evolution of the unstable type-I solution leads directly 
to the establishment of the stable type-II soliton, or, instead, transition to an 
excited state ( "breather" ) is observed, which then slowly relaxes into the stable 
type-II soliton. Also shown are regions where the stable soliton does not self- 
trap; instead, the unstable soliton either decays to zero, or generates a stable 
delocalized state. 

The solitons of types I and II are actually connected by a bifurcation. Indeed, 
the model which supports a stable soliton is a bistable system, as the zero 
solution is stable too in this case. According to well-known general principles, 
in a bistable system a separatrix must exist, which is a border between attraction 
basins of the two stable states. Usually, the role of the separatrix is played by 
an extra unstable soliton solution, whose amplitude is smaller and width larger 
than those of the stable soliton (see, e.g., Refs. 24 and |221)- The unstable 
type-I soliton is such a separatrix in the case oi q > qc = 0.25. 

The fact that the type-I soliton is stable at g < gc, while no stable soliton 
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Figure 8: The shape of the unstable type-I (separatrix) and stable type-II soli- 
tons in the FH (a) and SH (b) components for for a value of the mismatch 
parameter exceeding its critical value, q — 0.3 > Qc = 0.25. The other parame- 



ters are P — 0, ai = 0.3, and ao — {ao)e 



0.36, as per Eq. 
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Figure 9: The stability map for the type-II soliton self-trapping from the un- 
stable separatrix (type-I) soliton in the region q > 0.25. In the regions marked 
as "fast" and "slow" relaxation, the unstable soliton directly relaxes into the 
stable one, or does so via formation of a breather with slowly decaying intrinsic 
oscillations. In the region of "unstable solitons", the unstable soliton always 
decays to zero, and in the region of "delocalized solutions", a spatially periodic 
state sets in, instead of a stable soliton. The map pertains to /3 = 0, with ao 
chosen as per Eq. (|12|l . 
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of type II was found in that region, implies that a stabihty-swap bifurcation, 
involving both solitons, occurs at q — Qc- A conjectured bifurcation diagram is 
schematically depicted in Fig. ^| 




Stable solutions 



q ( detuning parameter ) 



Figure 10: The assumed stability-swap bifurcation between the solitons of types 
I and II (thin and bold curves) which occurs at the critical value of the mismatch, 
q ^ Qc = 0.25. Stable and unstable solutions are depicted by solid and dashed 
lines, and the vertical axis refers to the solitons' amplitudes. 

Note that the bifurcation diagram includes an unstable branch in the region 
of g < Qc, which is a conjectured continuation of the type-II solution to q < qc- 
This branch cannot be readily found from the stationary version of Eqs. JQ) and 
(|21l, since a good initial guess is not available for it, and it cannot be found in 
dynamical simulations either, as it is unstable. Thus, while we did not make a 
strong effort to explicitly find this branch - as, being unstable, it has no direct 
physical significance - we assume that such a branch exists. 

3.4 Stability of solitons in a generalized model 

The model based on Eqs. Q and Q can be generalized by replacing the 
evolution and spatial variables t and x, respectively, by z and r: 



This model corresponds, instead of the spatial solitons in cavities, to temporal 
solitons in waveguides. The temporal solitons are localized in the reduced-time 
variable r, and propagate along the coordinate z. In Eqs. (|14|) . the second 




(1 — iai) u + ianv* + il3 (u*) , 



(2q — iai) v + iaQU* . 



(14) 
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derivatives account for the temporal dispersion [rather than diffraction, in the 
original model Q, Q], D being a relative SH/FF dispersion coefficient. We 
have checked that, in this generalization (for instance, with D = 2, instead of 
D = 1), the results for the shape and stability of solitons are very similar to 
those reported above. 



3.5 Investigation of moving solitons 

A straightforward extension of the above analysis is to search for moving soli- 
tons. Their existence in the present model is not obvious, as the drive terms 
in Eqs. l|2Jl and 10} (unlike the loss terms) break the Galilean invariance of the 
equations. Without the drive, the Galilean boost with an arbitrary speed C, 



e 



iCX „, , „, 2iCX 



e^^^^, (15) 



transforms any solution u{x), v{x) into a moving one, u{x — Ct), v(x — Ct). 

To investigate the possibility of the existence of moving solitons in the 
present model, we ran systematic numerical experiments, taking the stable soli- 
ton solutions in the analytical or numerical form, as found above, and simulating 
the full system of Eqs. Q and ^ with the boosted initial conditions (|15|l . The 
result is that steadily moving solitons are not possible. Instead, a critical value 
Ccr of the boost parameter C was found, such that the soliton moves for a while 
but quickly comes to a halt and returns to its initial form, if C < Ccr, as shown 
in Fig. El In the opposite case, C > Ccr, the soliton always gets destroyed, see 
an example in Fig. 1)12(1 . 




Figure 11: An example of motion and subsequent stoppage of the boosted soliton 
®, in the case of ai = 0.02, q = 0.245, /? = 0, ao = (ao)cxact = 2-828 [see Eq. 
H12(l ]. and C = 0.3 (this value of the boost is close to the critical value, which 
is Ccr = 0.3335 in the present case). The panels (a) and (b) show the evolution 
of the pulse's amplitude profiles in the FF and SH components, respectively. 
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Figure 12: The same as in Fig. ^1 but with the boost parameter C = 0.35 
shghtly exceeding the critical value Ccr- 



The critical boost Ccr depends on values of the model's parameters. In 
particular, for the family of the exact solitons taken in the form of Eqs. 
with ao selected as per Eq. (|12|) . Cd- is shown, as a function of the loss parameter 
«! at different fixed values of the mismatch q, in Fig. (|13|l . 

The particular mode of the destruction of the boosted soliton in the case 
of C > Ccr depends on values of the parameters. In particular, instead of the 
straightforward decay, as in Fig. 1121 the soliton may first split into two secondary 
pulses with different velocities, both of which eventually decay. An example of 
the latter scenario is displayed in Fig. E| Also possible are situations in which 
the soliton does not split, but its FF component decays much faster than the 
SH one (although the FF and SH loss parameters are equal). 

4 Interactions between solitons 

Once stable solitons have been found, the next necessary step in the investigation 
of their fundamental dynamical properties, as well as in the development of 
potential applications, is the study of interactions between them. To this end, we 
ran systematic simulations of configurations initially composed of two identical 
stable solitons with centers placed at a distance Xq. Note that, unlike the 
situation in models of the Ginzburg-Landau type (the ones with the intrinsic 
gain) , in models with the parametric gain the phase of each individual soliton is 
locked to a single value [see Eqs. © and l|ll|l ]. hence the relative phase of the 
two solitons is not a free parameter, but is equal to zero 25 ; for this reason, 
the solitons always attract each other. The simulations demonstrate that, in all 
the cases, the attraction gives rise to merger of the two solitons into a single 
one. If the loss parameter ai is large enough, so that the soliton existing at 
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Figure 13: The critical value Ccr of the boost parameter (initial speed), which 
is defined in Eqs. H15|l . as a function of ai at fixed values of q, for the exact 
solitons (jHl [i.e., with /3 = 0, and = (Q!o)cxact]- ^he soliton comes to a halt, 
or gets destroyed, respectively, in the cases C < and C > Ccr, as shown in 
Figs. CHandini 




X 



(a) (b) 

Figure 14: An example of the splitting of the soliton ||HJ) to which an overcritical 
boost was applied, and subsequent decay of the secondary pulses. In this case, 
ai = 0.15, q = 0.2, /3 = 0, ao = (ao)oxat = 0.25, the critical value of the boost 
is Ccr ~ 0.7, and the actual value of the boost is C = 0.75. The panels (a) and 
(b) show the evolution of the amplitude profiles in the FF and SH components, 
respectively. 
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this value of ai has no complex eigenvalue of small intrinsic perturbations (see 
Figs. 0]and|5l), the resulting single soliton emerges in its stationary form. On 
the other hand, if ai is small, and the soliton possesses an intrinsic complex 
eigenvalue, the final soliton appears in an excited (vibrating) state, which then 
slowly relaxes to the static one. The time necessary for the fusion of the two 
solitons into one depends on the initial separation Xq, but the outcome of the 
interaction does not depend on Xq. 

Below, we illustrate these conclusions by typical examples. In all the cases, 
we used exact soliton solutions (jSJ to construct the initial configuration. Sim- 
ulations of the configuration composed of two solitons of a more general form 
produced virtually the same results as those displayed below. 

If the loss parameter ai is small enough, a pair of exact solitons merge into 
a pulse which demonstrates nearly persistent intrinsic vibrations, as it was said 
above and is shown in Fig. 1151 Continuing the simulations on a much longer time 
scale demonstrates that the vibrations slowly fade out, and the pulse relaxes to 
the static configuration. The latter feature is seen clearer in Fig. ^1 which 
corresponds to a smaller initial separation between the solitons. 




Figure 15: The merger of two solitons, with the initial separation Xq = 5.2 
between them, into a pulse with excited internal vibrations, which later slowly 
relaxes into the stationary soliton. The panels (a) and (b) show the FF and SH 
components of the field, respectively. The parameters are q = 0.207, ai = 0.094, 
and ao = (ao)exact = 0-196, see Eq. 

If ai is larger, so that the soliton does not support intrinsic oscillatory modes, 
two solitons, even separated by a relatively large distance, fuse into a single 
soliton which emerges in the stationary state (without intrinsic vibrations), as 
shown in Fig. El The fact that the final soliton is identical to each initial one 
is obvious from Fig. ^1 which compares the initial field configuration in both 
component and its eventual shape. 
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Figure 16: The same as in Fig. ^1 but for a smaller initial separation between 
the solitons, Xq — 3.6. In this case, it is obvious that the vibrating pulse, 
produced by the merger of the two initial solitons, relaxes towards the static 
soliton. 
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(a) (b) 

Figure 18: The comparison of the initial and final wave-field profiles in the FF 
(a) and SH (b) components in the same case which is shown in Fig. El 

5 Conclusion 

In this work, we have introduced a model of a lossy second-harmonic-generating 
(t^(2)^ cavity driven by a pump wave at the third harmonic, which gives rise 
to a new type of driving terms, characterized by the cross-parametric gain. 
The equation for the fundamental-frequency wave may also contain a quadratic 
self-driving term, which is generated by the x^'^^ nonlinearity. 

Unlike previously studied phase-matched models of x*-^^ cavities driven through 
down- or upconversion, the present model admits the exact analytical solution 
for the soliton, at the specially chosen value of the gain parameter. Two general 
families of soliton solutions were found in a numerical form, one of which is a 
continuation of the exact analytical solution. At given values of the parameters, 
one soliton is stable, and one is not. They swap the stability via a bifurcation 
which occurs at a critical value of the mismatch parameter. Full stability regions 
of the solitons were identified by means of numerical computation of the corre- 
sponding eigenvalues for small perturbations (stability conditions for the zero 
solution, which is a necessary ingredient of the full conditions for the soliton's 
stability, were found in an analytical form). The stability of the solitons was 
also verified in direct simulations, with the conclusion that the unstable soliton 
rearranges into the stable one (which may appear in the form of a breather), 
or into a delocalized state, or decays to zero. Additionally, it was found that 
steadily moving solitons do not exist in the present model. If the soliton is ini- 
tially boosted, it either comes to a halt, or, if pushed too hard, gets destroyed 
(possibly, via splitting into two pulses). 

Interactions between initially separated solitons were also investigated by 
dint of systematic direct simulations. It was found that stable solitons always 
merge into a single one. In the system with weak loss, the final solitons ap- 
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pears in an excited form (the breather), and then slowly relaxes to the static 
configuration. If the loss is stronger, the final soliton emerges in the stationary 
form. 

The model introduced in this work can be further investigated in various 
directions. In particular, a two-dimensional version of this cavity model may be 
an interesting subject. 
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